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Symbolic dynamics and chaos in plane Couette flow 

Y. Charles Li 


Abstract. According to a recent theory [5], when the Reynolds number is 
large, fully developed turbulence is caused by short term unpredictability 
(rough dependence upon initial data); when the Reynolds number is mod¬ 
erate, often transient turbulence is caused by chaos (long term unpredictabil¬ 
ity). This article aims at studying chaos in plane Couette flow at moderate 
Reynolds number. Based upon the work of L. van Veen and G. Kawahara [8] 
on a transversal homoclinic orbit asymptotic to a limit cycle in plane Cou¬ 
ette flow, we explore symbolic dynamics and chaos near the homoclinic orbit. 
Mathematical analysis shows that there is a collection of orbits in the neigh¬ 
borhood of the homoclinic orbit, which is in one-to-one correspondence with 
the collection of binary sequences. The Bernoulli shift on the binary sequences 
corresponds to a chaotic dynamics of a properly defined return map. 


1. Introduction 

From our analytical result |5j (supported by our numerical experiments [3]), 
perturbations in turbulence amplify in time according to e CT vfle\/t+o-it w ]- lere (Jl _ 
and a depends only on the base solutions on which the perturbations are 
introduced. When the time is small, the first term in the exponent dominates, and 
this term can cause the amplification to be superfast when the Reynolds number 
is large. By the time t ~ Re, the two terms in the exponent are about equal. 
After the time t ~ Re , the second term dominates, and this term is the classical 
Liapunov exponent that causes chaos (long term unpredictability). Thus the time 
t ~ Re is the temporal separation point between short term unpredictability and 
long term unpredictability. When the Reynolds number is large, long before the 
separation point t ~ Re, the first term in the exponent already amplifies the per¬ 
turbation to substantial size, and the second term does not get a chance to act. 
Thus fully developed turbulence is dominated by such short term unpredictability. 
When the Reynolds number is moderate, both terms in the exponent have a chance 
to dominate, and the corresponding (often) transient turbulence is dominated by 
chaos in long term. Since random perturbations can happen at any time, this leads 
to the fact that the exact mathematical solution of the Navier-Stokes equations 
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does not describe any turbulent flow even approximately in short time, when the 
Reynolds number is large. This is in sharp contrast to chaos for which the exact 
mathematical solution at least approximates a chaotic flow in finite time. We study 
the large Reynolds number Navier-Stokes equations somewhere else [8]. Here we 
study moderate Reynolds number Navier-Stokes equations and focus upon sym¬ 
bolic dynamics and chaos near the homoclinic orbit discovered by L. van Veen and 
G. Kawahara [§]. There have been a lot of dynamical system studies on moderate 
Reynolds number Navier-Stokes equations, and various dynamical system objects 
have been discovered pQ, e.g. fixed points, periodic orbits, periodic doubling bi¬ 
furcation heteroclinic tangle [7j. These studies were conducted on specific 

flows like plane Couette flow, plane Poiseuille flow, and pipe Poiseuille flow. 


2. Bernoulli shift on binary sequences 

The Bernoulli shift is an automorphism that has chaotic dynamics. In fact, it is 
a canonical chaos map to which other chaos maps can be shown to be topologically 
conjugate to. Let E be the set of all doubly infinite binary sequences: 

a = (• • • a_ 2 a_ia 0 , aia 2 • • •), 

where a n = 0 or 1 for any integer n. For any point in E, 

a = (■■■a_ 2 a_ 1 a 0 ,a 1 a 2 ---), 
one can define the neighborhood basis 

Uj(a*) = {a € E | a n = a*, \n\ < j} 

where j = 1, 2, - - -. These Uj(a*)’ s generate a topology for E. Endowed with 
this topology, E becomes a topological space. The larger the j is, the smaller the 
neighborhood Uj(a*) is. The Bernoulli shift x is a map defined in E. For any point 
a in E, let b = x(a) , then b n = a n +i for all integer n. That is, if a is given by 

a — (• • • a_ 2 a_iao, aia 2 • • •), 

then x( a ) is given by 

y(a) = (• • • a_ 2 a_ia 0 ai,a 2 • • •), 

in other words, the comma is moved one step forward. The inverse Bernoulli shift 
is simply moving the comma one step backward. Both the Bernoulli shift and its 
inverse are homeomorpisms. Thus, the Bernoulli shift is an automorphism. The 
Bernoulli shift has sensitive dependence upon initial data - the signature of chaos. 
Take two close initial points a and b in E, say b is in the neighborhood Un{o) for 
a large TV, but not in f7/v + i(a), that is 

a n = b n , \n\ < N, 


but 


either oat ^ bpj or a-w ^ fo-vr- 


After applying either the Bernoulli shift or its inverse N times, the images are not 
inside Ui neighborhood, that is, either x N (b) is not in Ui(x N ( a )) or X~ N (b) is not 
in U\{x~ N (a))- Thus no matter how close two different initial points are, after 
enough iterations of the Bernoulli shift, the images are far away. 
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3. van Veen - Kawahara transversal homoclinic orbit in plane Couette 

flow 

For the plane Couette flow at moderate Reynolds number (~ 400), the topo¬ 
logical setup of the homoclinic orbit numerically observed by L. van Veen and G. 
Kawahara .8] can be depicted as in Figure [I] where H is the homoclinic orbit as¬ 
ymptotic, in both forward and backward time, to the limit cycle (isolated periodic 
orbit) P. The limit cycle P has a two-dimensional unstable manifold and a one- 
codimensional stable manifold, which intersect transversally along the homoclinic 
orbit H. In [5], the limit cycle is very small, and the homoclinic orbit is very large. 
Nevertheless, the topology is the same as in Figure [I] In both forward and back¬ 
ward time, the homoclinic orbit just keeps wrapping around the limit cycle closer 
and closer. From the analysis of partial differential equations, the Poincare return 
map is not well-defined on the Poincare section transversal to the limit cycle. Thus 
one cannot reduce the problem of establishing chaos near the homoclinic orbit to 
that on the Poincare section, instead we establish chaos directly in the original 
space [4j. 



Figure 1. Topological setup of the homoclinic orbit. 


4. Pseudo-orbits 

To establish the symbolic dynamics and chaos near the homoclinic orbit, we 
introduce pseudo-orbits |4j. Let c be a point on the limit cycle P, and ryo be the 
periodic orbit starting from c, that travels on P many times (i.e. many copies of 
P). We can follow the homoclinic orbit backward until it wraps around the limit 
cycle many times and reach a point a that is close to the point c on the limit cycle. 
Then the point a is very close to the limit cycle. Also we can follow the homoclinic 
orbit forward until it wraps around the limit cycle many times and reach a point 
b that is close to the point c on the limit cycle. Again the point b is very close to 
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the limit cycle. See Figure [2j To make the illustration in Figure [2] clearer, we did 
not draw the wrapping of the homoclinic orbit around the limit cycle. From b to 
a, we put a mollifying connector that is tangent to the homoclinic orbit at b and a, 
and tangent to the limit cycle at c. With this connector on, we get a closed loop 
denoted by ift. While rji is not an orbit, it is a pseudo-orbit since the jump from b 
to a is small. The starting point and the end point of the pseudo-orbit i]± is c. To 



Figure 2. Pseudo-orbit setup. 

any doubly infinite binary sequence 

a = (■ ■ ■ a,- 2 a -i a ch a i °2 • • • )> 

where a n = 0 or 1 for all integers n, there corresponds a pseudo-orbit 

V a = ’ Vo,— 2 Vo.— i Vao ) Va 1 Va 2 ’ ' ’ )• 

Note that the connecting point between ijo and 771 is c. 

5. Symbolic dynamics and chaos 

Through contraction map argument on the graphs of pseudo-orbits [41, we can 
show that to any doubly infinite binary sequence 

a = (• • • a_ 2 a_ia 0 , aia 2 • • •), 

where a n = 0 or 1 for all integers n, there corresponds an true orbit O a that is close 
to the pseudo-orbit r] a [3J (strictly speaking, the analysis in |4j applies to all the 
Galerkin truncations of the Navier-Stokes equations, which are the systems that 
numerical simulations truly calculate. Due to analytical technicality, the analysis 
in 3] does not directly apply to the full Navier-Stokes equations due to the singular 
perturbation nature of the Laplacian dissipative term). We call O a the shadowing 
orbit of the pseudo-orbit rj a . Thus the topological space E of all doubly infinite 
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Figure 3. Poincare section setup. 



Figure 4. Cantor subset of the Poincare section. 

binary sequences is in one-to-one correspondence with the set S of these true or¬ 
bits Oa s. We put a Poincare section S in a location before the homoclinic orbit 
starts to wrap around the limit cycle in backward time (see Figure [3]). Let h be 
the intersection point of homoclinic orbit with the Poincare section S before the 
homoclinic orbit starts to wrap around the limit cycle in backward time. Let p be 
the intersection point of limit cycle with the Poincare section S. h and p are far 
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away apart. See Figure |3J We designate h the chosen intersection point of the loop 
with the Poincare section S 7 and p the chosen intersection point of the first copy 
of the limit cycle P in the loop 770 . We designate the chosen intersection point of 
p a with the Poincare section S to be that of p ao - Finally, we designate the chosen 
intersection point of O a with the Poincare section S to be the one that is close to 
that of p a . See Figure [4] for an illustration. Denote by A the set of intersection 
points of all O a ’s with the Poincare section S for all doubly infinite binary sequences 
a £ E. A has a Cantor type structure. The topology of A discussed above induces a 
topology on A. For any point o* £ A which is the intersection point of O a . with the 
Poincare section S, the neighborhood I7j(a*) of a* induces a neighborhood W)(o*) 
of o*, which consists of the intersection points of all O a with the Poincare section 
S for all a £ The Poincare return map F is not defined on the Poincare 

section S as an analytical result of partial differential equations, but it is defined on 
A. Defined on the topological space A, the Poincare return map F is topologically 
conjugate (equivalent) to the Bernoulli shift y defined on the topological space E. 
That is, there is a homeomorphism <f> : E —> A such that the following diagram 
commutes: , 


E 


A 


*1 i F 

E —» A 

Since the Bernoulli shift y has sensitive dependence upon initial data in the topo¬ 
logical space E, the Poincare return map F also has sensitive dependence upon 
initial data in the topological space A, that is, the dynamics of F is chaotic. If 
we choose two different nearby points in A, then after enough iterations or inverse 
iterations of F, the image of one point is close to p and that of the other is close to 
h (note that p and h are far away apart). 
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